Coexistence of coupled magnetic phases in epitaxial TbMnOs films revealed by 

ultrafast optical spectroscopy 
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Ultrafast optical pump-probe spectroscopy is used to reveal the coexistence of coupled antiferro- 
magnetic/ferroelectric and ferromagnetic orders in multiferroic TbMnOs films through their time 
domain signatures. Our observations are explained by a theoretical model describing the coupling 
between reservoirs with different magnetic properties. These results can guide researchers in cre- 
ating unique kinds of multiferroic materials that combine coupled ferromagnetic, antiferromagnetic 
and ferroelectric properties in one compound. 



Multiferroic materials have attracted much interest in 
the past decade, due not only to their potential device ap- 
plications, but also their manifestations of coupling and 
interactions between different order parameters (particu- 
larly electric polarization and magnetic order). Recently, 
much attention has been focused on pcrovskite mangan- 
ites, i?MnC>3 (i?=rare earth ions), due to the discovery 
of a large magnetoelectric effect in these materials [H, [J. 
The first member of this family to be discovered was 
TbMnC>3 (TMO), which is now well established as a typ- 
ical magnetoelectric multiferroic. 

Extensive studies have already been done on single 
crystal TMO (SC-TMO) [HQ. In brief, SC-TMO, with a 
distorted orthorhombic perovskite structure, has an anti- 
ferromagnetic (AFM) phase transition at Tn — 40 K with 
sinusoidally ordered Mn moments. Below Tfe — 28 K, 
ferroelectric (FE) order develops owing to the appearance 
of cycloidal spiral AFM spin structure. In contrast, there 
are relatively few reports describing the properties of 
TMO thin films (typically grown on SrTi0 3 (STO) sub- 
strates) [H-IH . In general, thin films can enable enhanced 
functionality in materials, as their physical parameters 
can be changed by modifying their structure via strain 
imposed by the substrate. Strain, in particular, can pos- 
sibly provide a different route to direct FE-FM coupling, 
which is very rare @, Q . This could benefit electronic de- 
vice applications by providing low power consumption, 
high speed operation, and greater electric/magnetic field 
controllability |{|. 

In TMO films, previous investigations of magnetic 
properties revealed an unexpected ferromagnetic (FM) 
order [1,0,0], in contrast to SC-TMO. However, none 
of these reports have demonstrated the existence of FE 
(or equivalently, spiral AFM) order in these thin films. 
Furthermore, one can ask: if FE order does exist, is it 
coupled to FM order? 

Here, we use ultrafast optical pump-probe spec- 
troscopy to reveal coexisting coupled magnetic and fer- 
roelectric orders in epitaxially grown TMO thin films, 
which were not observed in previous work [3 8] . Our tem- 



perature ( T)-dependent transient differential reflectivity 
(AR/R) measurements show clear signatures of FE and 
AFM phases developing as the film thickness changes. 
We propose a model to explain the interactions between 
AFM/FE and FM orders. These results reveal that the 
coupling between different magnetic and FE orders ob- 
served in our TMO thin films may offer greater control 
of functionality as compared to bulk multiferroics. 

Our transient differential reflectivity experiments used 
a Ti: Sapphire laser oscillator with a repetition rate of 80 
MHz producing ^140 femtosecond (fs) pulses at the sam- 
ple position, with a center wavelength of ^830 nm (1.5 
eV). The pump and probe pulses are cross-polarized and 
normally incident on the sample surface. The pump has 
a typical fiuence of 1 /iJ/cm 2 , which is 10 times larger 
than that of the probe. This excites ~ 10 -4 quasiparti- 
cles/unit cell on average, ensuring that the photoexcita- 
tion does not induce large thermal or non-thermal effects. 

The TMO films investigated here were grown on 
STO(001) substrates using pulsed laser deposition 
(PLD), with thicknesses of <i=34, 50, and 150 nm. 
We also examined a TMO single crystal (SC-TMO) 
grown through the floating zone technique and a 50 nm 
SmMn0 3 (SMO) film grown by PLD on a STO (001) 
substrate. In the bulk, SMO displays A-type AFM order 
below ~ 60 K, but unlike TMO, does not exhibit sinu- 
soidal/cycloidal spin order or FE order [l(J. Therefore, 
comparing measurements on our TMO films with those in 
a TMO single crystal and a SMO film will help us unravel 
the influence of these order parameters on quasiparticle 
dynamics in TMO. 

Figures [Ha), (b) and (c) show typical AR/R data 
obtained on a d =50 nm TMO film. Very similar re- 
sults are obtained on the 34 and 150 nm films, except 
that the initial signal magnitude increases with thick- 
ness (Fig. [U(cl)). The dynamics measured in the SMO 
film (Figs. QTe) and (f)) are also quite similar to those 
measured in the TMO films. We observe three distinct 
relaxation processes in these films, occurring on different 
timcscales. Within a time t ~1.5 picoseconds (ps), there 
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FIG. 1. AR/R signals measured on a 50 nm TMO film on 
long ((a) and (c)) and short timescales (b). Part (d) shows 
ultrafast dynamics at 20 K on short timescales in three TMO 
films with d ~34 nm (dashed), 50 nm (solid), and 150 nm 
(dashed-dotted) , respectively. The black solid curve in (a) 
demonstrates the fitting. For comparison, data at 20 K for 
SC-TMO (blue) is also shown in (c) and (d). AR/R signals 
measured on a 50 nm SMO film on long and short timescales 
are shown in (e) and (f), respectively. 



is an ultrafast increase in the magnitude of the reflec- 
tivity, followed by a fast signal reversal. The former is 
due to the photoexcitation of carriers via intersite Mn 
d-d transitions, which dominate the optical response at 
1.5 eV in perovskite manganites [111, |l2|. We note that 
the strong temperature dependence observed in our data 
further supports this, and indicates that p-d transitions 
do not play a significant role [l2j]. The fast signal rever- 
sal, occurring at t ~ 1 ps, can be mainly attributed to 
electron-phonon (e-ph) relaxation, as in other mangan- 
ites [IH . The timescale of these initial ultrafast changes is 
approximately independent of the temperature and film 
thickness. 



After the initial ultrafast dynamics are complete, an 
additional rising process (in which the magnitude of 
AR/R increases, indicated by the arrows in Fig. [TJ, 
clearly appears in our films around ^1.5 ps below a tem- 
perature Tc (~140 K) followed by a much slower relax- 
ation process, with a decay time of hundreds of picosec- 
onds, that returns the photoexcited system back to equi- 
librium. The latter process is likely due to heat diffusion, 
as will be discussed in more detail below. The rising 
process, occurring on an intermediate time scale of tens 
of ps at low temperatures (T), is unlikely to be due to 
long-range bulk-like AFM or FE order (as one would ex- 
pect in SC-TMO) for two reasons. First, the transition 



temperature for the emergence of this additional rising 
process, T c , is well above T N and T FE of SC-TMO (Fig. 
QJb)). Secondly, our data reveals that this process is ab- 
sent in SC-TMO (Figs. QTc) and (d)). Instead, we only 
found two relaxation processes in SC-TMO, which have 
timescales of ^1 ps and hundreds of picoseconds, respec- 
tively. The short-lived ps dynamics in SC-TMO, repre- 
senting rapid thermalization of the excited charge, lattice 
and spin subsystems, are consistent with previous mea- 
surements on other long-range AFM ordered compounds 
l~3l- 16| . an d the slower dynamics arc due to thermal dif- 
fusion !(14f . The rising process observed in our TMO (and 
SMO) films is thus unlikely to originate solely from the 
properties of their bulk counterparts. Instead, we note 
that a relatively slow spin-lattice relaxation process asso- 
ciated with long-range FM order was previously observed 
in some perovskite oxides [liiEIIj- Therefore, the rising 
process observed here might be linked to the presence of 
FM order in our films. 
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FIG. 2. (a) The magnetic moment versus magnetic field H at 
T = 5 K is shown for a 50 nm TMO film. ZFC-FC curves are 
shown in (b). The inverse susceptibility (x" 1 ) °f the nml is 
plotted as a function of temperature, T, in (c). The data was 
taken with a 1 kOe field applied parallel to the sample plane. 

To investigate this possibility, we measured the field 
(TJ)-dependent magnetic moment in our films using a 
SQUID magnetometer. Figure HJa) shows hysteretic be- 
haviour in the 50 nm TMO film, confirming the presence 
of ferromagnetism. This is supported by temperature- 
dependent measurements (Fig. Efb)), where a clear 
splitting of field-cooled (FC) and zero-field-cooled (ZFC) 
curves below T* ~48 K reveals the appearance of a net 
magnetization. To uncover more information on the mag- 
netic order (s), we plotted the inverse susceptibility x" 1 
as a function of temperature (Fig. Hfc)). We did not 
observe clear anomalies at x _1 (Tf,e) m our TMO films. 
By fitting the data in the high-temperature regime with 
the Curie- Weiss form x _1 oc (T — 6cw), we found the 
Curie- Weiss temperature, 6cw — —140 K. The nega- 
tive Ocw indicates that the exchange interactions in our 
films are primarily antiferromagnetic flij . Similar trends 
were observed in our other films, consistent with previ- 
ous results 0, S, S|- However, we note that the exact 
mechanism causing the observed ferromagnetism is still 
unclear. For example, ferromagnetism could arise from 
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a strain-induced change in the balance between the dif- 
ferent magnetic exchange interactions 0, H|, electronic 
phase separation (due to nucleation of a ferromagnetic 
phase that is induced by epitaxial strain pol]). or the 
presence of twin domain walls due to film strain relax- 
ation 
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Our magnetization measurements thus indicate that 
the rising process is linked to the presence of FM order in 
our films; however, this process emerges at Tc, well above 
T*, suggesting that additional physics could be involved. 
For more insight, we fitted our data for t > 1 ps with 
the equation: A r e~ t/Tr + A d e~ t/Td (A r > and A d < 0) 
(Fig. QJa)). This revealed some surprising anomalies in 
the temperature variation of the intermediate (r r ) and 
slow (ry) time constants for the SMO and TMO films 
studied here (Figure [3]) that were not observed in the 
SQUID measurements. 

In particular, very strong peaks in both r r (T) and 
r d (T) are found in the SMO film at -60 K. As dis- 
cussed above, SC-SMO has an A-type magnetic struc- 
ture below the Neel temperature (T/v — 60 K). There- 
fore, the anomalies in the measured time constants at 
—60 K should be associated with the appearance of long- 
range bulk-like AFM order. Similarly, in our TMO films, 
both r r and T d show clear kinks at —30 K and —46 K 
for e?=50 nm and 150 nm (most obvious in the latter), 
as indicated by the arrows. If either of these anomalies 
is primarily due to FM order, we also expect to observe 
it in the 34 nm TMO film based on the SQUID results 
and strain-related ferromagnetism. However, no obvious 
anomalies were found for the 34 nm film. Thus, as in 
SMO, since 30 K and 46 K are very close to Tfe and 2~W 
in SC-TMO, the anomalies in r r and r d can be associ- 
ated with the bulk-like FE and AFM phase transitions, 
respectively. In addition, the fact that 46 K is so close 
to T* suggests that TMO films may undergo AFM and 
FM transitions at the same temperature. 

It is relatively straightforward to explain the anomalies 
in the slower time constant, r d . This process originates 
from thermal diffusion, including both lateral diffusion 
within the films as well as diffusion into the substrate 
(the former becomes more important with increasing film 
thickness). This can be shown by considering the diffu- 
sion equation dT/dt = aV 2 T, in which the heat transfer 
rate, reflected in l/r d , decreases with the thermal diffu- 
sivity a (= n/Cp), where k, C and p are the heat con- 
ductivity, the specific heat and the density, respectively. 
Thus, when and C peak at the magnetic phase tran- 



sition temperatures for TMO and SMO (l0j,|22|, thermal 
diffusion will correspondingly slow down, explaining the 
anomalous peaks in r d . 

The origin of the intermediate time constant r r re- 
quires more detailed consideration. In all TMO films 
studied here, we found Tc — \0cw\, suggesting that 
not only FM ordering, but also AFM exchange inter- 
actions, may influence this rising process. This could 
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FIG. 3. Decay time constants (r r and Td) at different temper- 
atures for the d=50 nm SMO film ((a) and (b)), and for d=34 
nm, 50 nm, and 150 nm TMO films ((c) and (d)). For clarity, 
TV for the 50 nm and 150 nm TMO films in (c) is shifted up by 
3 ps and 6 ps, respectively. Td for the 50 nm and 150 nm TMO 
films in (d) is shifted up by 50 ps and 150 ps, respectively. 
The arrows indicate the phase transition temperatures. 



occur through strong magnetic frustration due to com- 
peting magnetic interactions, which is often observed in 
multiferroics 0. One signature of this is when \0cw\ 
is much greater than T* and Tjy [H, [lij], which is the 
case for our films (Fig. [2]). This is further supported by 
our ultrafast optical data, where clear anomalies associ- 
ated with other magnetic orders are found in t t at T/v 
and Tfe, indicating that they coexist with the FM order 
revealed by our SQUID measurements. The magnetic 
interactions inducing the long-range bulk- like AFM/FE 
order may thus compete with the magnetic interactions 
contributing to ferromagnetism, influencing quasiparticle 
dynamics on the timcscale of r r . 

Motivated by these findings, we propose a model with 
two spatially separated magnetic reservoirs, one with 
bulk- like sinusoidal/spiral AFM order and the other with 
FM order, to describe the intermediate rising dynamics 
observed in i?MnO,3 films. They are characterized by two 
effective temperatures, T A and T F , respectively. The 
coupling between the two reservoirs occurs via the ex- 
change spin-spin interaction Ja-fS a ■ S F between these 
two systems [23j . This coupling in the photoexcited sys- 
tem is manifested through energy transfer via magnon- 
magnon coupling/conversion between the two reservoirs, 
where the magnon excitations in both reservoirs are in- 
troduced through intersitc d-d transitions [HI]. Similar 
to the classic two-temperature model describing electron- 
phonon coupling 24| , the dynamical behavior of the pho- 
toexcited coupled system can be described by the time 
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evolution of its effective temperatures, given by 

C A j t {T A ) =g A - F (T F -T A ), 

C F j t (T F )=9A- F (T A -T F ), (1) 

where §a-f is the energy transfer rate, and C J (j = 
A, F) are the magnetic heat capacities for the two reser- 
voirs with AFM order and FM order, respectively. After 
solving these differential equations, we found that the 
non-equilibrium coupled system can be characterized by 
a time constant r, 

T ~' = 9A-F (j^a+^f)- ( 2 ) 

Now, based on this theoretical r that describes the in- 
termediate rising process, we can understand how the 
magnetic phase transitions could be revealed in the tem- 
perature dependence of r r (Fig. |3|). Specifically, r will 
reach a maximum in the vicinity of a critical tempera- 
ture because: (1) the heat capacities C A and C F peak at 
critical temperatures (Io| . and (2) qa-f might become 
smaller at the magnetic phase transition temperatures 
due to the critical slowing down phenomenon [251 ] . Fur- 
thermore, if we assume C A and C F have similar values, 
we can estimate the coupling constant gA-F — 2 x 10 10 
WK^mole- 1 for TMO films at 30 K, calculated by tak- 
ing the spin specific heat C A ~ 1 JK _1 molc _1 and 
r r ~ 25 ps [HI [2^. Our model thus shows that the 
anomalous peaks in r r {T) are due to the presence of co- 
existing, coupled magnetic orders in these films. 
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FIG. 4. Temperature dependence of A" 1 for three TMO films. 

It follows that the appearance of the intermediate ris- 
ing process (A r > 0) below \9cw \ represents the competi- 
tion between magnetic interactions. Thus, the amplitude 
A r reflects the coupling between neighboring spins, ex- 
pressed by the spin-spin correlation function < SiSj >, 
which includes intra-reservoir correlations between neigh- 
boring spins in each reservoir as well as the inter-reservoir 
correlation function < S A S F >; this last term also con- 
tributes to the coupling constant gA-F in Eq. (fT]). Fig- 
ure Q] demonstrates that A~ 1 decreases with decreasing 



T. This may indicate that < SiSj > increases as the 
temperature decreases, as in previous work [ll|. Most 
interestingly, we found a very clear transition at T ~ 30 
K, comparable to the bulk Tfe, for both 50 nm and 
150 nm TMO films. Below ~ 30 K, A' 1 is nearly con- 
stant with temperature, unlike its behavior for T > 30 
K. Since A r (T) is influenced by the spin-spin correlation 
function < SiSj >, this result is consistent with previ- 
ous work showing that the propagation spin wave vector 
k s is locked at a constant value below Tfe (^30 K) [l[ 
when TMO undergoes a FE phase transition and the spin 
structure changes from sinusoidal to spiral; this further 
supports the existence of FE order in our TMO films. 

In conclusion, we have used ultrafast optical pump- 
probe spectroscopy to reveal coexisting, coupled mag- 
netic and ferroelectric orders in epitaxial TMO films 
through their time domain signatures. This could lead to 
an enhanced magnetoelectric effect in these films as com- 
pared to their bulk counterparts, which would allow the 
FE polarization and/or magnetization in TMO films to 
be more easily controlled by external fields than in bulk 
TMO. Our results therefore indicate that thin perovskitc 
manganite films could be an unique kind of multiferroic 
combining coupled FM, AFM and FE properties in a sin- 
gle material. 
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